Proposal for a geophysical search for dilatonic waves 
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Abstract 

We propose a new method of searching for the composition-dependent dilatonic waves, predicted 
by unified theories of strings. In this method, Earth's surface-gravity changes due to translational 
motions of its inner core, excited by dilatonic waves, are searched for by using superconducting 
gravimeters. This method has its best sensitivity at the frequency of ~ 7 x 10~ 5 Hz, which is lower 
than the sensitive frequencies of previous proposals using gravitational-wave detectors: ~ 10 to 
1000 Hz. Using available results of surface-gravity measurements with superconducting gravimeters 
and assuming a simple Earth model, we present preliminary upper limits on the energy density of a 
stochastic background of massless dilatons at the low frequency. Though the results are currently 
limited by the uncertainty in the Earth model, this method has a potential of detecting dilatonic 
waves in a new window. 
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I. INTRODUCTION 



The scalar gravitational fields (or the dilatons) appear in scalar-tensor theories of grav- 
itation and also in theories toward the unification of fundamental forces in nature, such 
as Kaluza-Klein, supergravity and string theories. Scalar-tensor theories must agree with 
general relativity within an accuracy of ~ 10~ 4 - 10~ 5 in the post-Newtonian limit, because 
of the results from the Cassini time-delay experiment [11 and the search for the Nordtvedt 

n n 

effect (i.e. violations of the universality of free-fall for massive bodies [2[]) [3|. However, 
in strong gravity systems, in which the post-Newtonian approximation is not applicable, 
they could exhibit large observable deviations from general relativity |4( and other tests 
involving a strong gravity source are necessary. One such tests is to search for scalar waves 
from a spherical symmetric gravitational collapse of a star 5), which are predicted by 
scalar-tensor theories but not by general relativity. 

In many models of unified theories, the scalar field mediates a new force (or a fifth force 

n 

[7J]), which leads to a violation of the universality of free-fall and/or a deviation from the 
inverse-square law of gravitation. Such deviations from the inverse-square law and violations 
of the universality of free-fall have been searched for and constrained by various experiments 
at different ranges (see e.g. ?], 0, However, string theory, one of the most promising 
theories for the unification, predicts the existence of a relic background of the dilaton, which 
could be a significant component of the dark matter in the present universe [l[J ; though the 
fifth-force searches indicate that the dilaton coupling to matter is very weak, the weakness 
of the coupling could be compensated by the high intensity of the background. Therefore, 
the detection of the scalar waves or a direct experimental constraint on the amplitude of the 
scalar waves would play an important role for the development of the unified theories and 
verification of theory of gravitation. 



Motivated 
( see e.g. jolfll 



jy these predictions, strategies of detecting the scalar waves have been studied 
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191 ] ) . These analyses show that some of the predicted 



scalar waves could be detected by using currently operating or advanced planned ground- 
based gravitational wave detectors. However, no experimental results have been presented 
yet. Here, we propose a new approach of searching for scalar waves, using superconducting 
gravimeters, and we tentatively estimate upper limits on the energy density of a stochastic 
background of dilatonic waves. The concept of this method (described in detail later) is 
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similar to the one of the geophysical test of the universality of free-fall [20] in the way 
that both methods use the Earth as the test body and superconducting gravimeters as the 
displacement sensor. 

Dilatonic waves can interact with detectors in two ways [2l|]: (i) directly, through the 
effective dilatonic charges of the detectors, which depend on the internal composition of the 
detectors, and (ii) indirectly, through the geodesic coupling of the detectors to the scalar 
component of the metric fluctuations. In the former case, the response of the detectors to 
the dilatonic waves is nongeodesic [21]. The proposed method is to seek for an effect of the 
direct coupling. 

The sensitive frequency range of ground-based gravitational wave detectors is about 10 
to 1000 Hz, while it is about 7 x 10~ 5 Hz for the proposed method; a new window can be 
explored with this method. According to a cosmological scenario motivated by string theory, 
a relic background of light dilatons, with mass as small as ~ 10~ 19 eV or less (depending 
on the level of the flatness of the background spectrum), could be a significant component 
of the dark matter (see Section 6 of 10] and references therein). The sensitive frequency 
range of the proposed method corresponds to the mass range of ~ 3 x 10~ 19 eV; the possible 
cosmological scenario could be tested with this method. However, in this paper, we focus 
on the massless dilaton for simplicity. 



II. THE CONCEPT 



In this method, we consider the Earth's inner core as the receiver of dilatonic waves. 
The inner core, enclosed in its liquid outer core, is weakly coupled to the rest part of the 
Earth mainly by gravitational forces. When a gravitational wave impinges on the Sun and 
the Earth, from any direction that is not aligned with the Sun-Earth line, their proper 
separation oscillates. Gravitational waves couple to matter in the universal way; to first 
order, the propagation of gravitational waves would not cause relative motions between the 
inner core and the rest part of the Earth. However, when dilatonic waves pass, because of the 
difference in dilatonic charge (namely, the difference in their chemical composition), there 
would be relative motions between the inner core and the rest part of the Earth (see next 
section for detail). Such relative motions would result in changes of the gravitational field at 
the surface of the Earth. The surface-gravity changes can be searched for by superconducting 



3 



gravimeters, which are the most sensitive instruments for measurements of gravity changes 
at low frequencies. 

The dilatonic charges of the inner core and the rest part of the Earth are different because 
of the difference in their chemical composition (to be discussed in detail later); the inner 
core mainly consists of iron and nickel, while the rest part of the Earth is mainly made of 
lighter elements such as silicon oxides. 

The response of the inner core to the dilatonic waves could be greatly enhanced when 
the waves have the same frequencies as the natural oscillation frequencies of the inner core. 



III. RESPONSE OF THE INNER CORE TO DILATONIC WAVES 

For simplicity, we assume a simple Earth model and configuration as shown in Fig. [IJ A 
dilatonic wave propagates along the Sun-Earth line, which connects the centers of mass of 
the Sun and the Earth. Their rest separation is L « 1.5 x 10 11 m. 

The Earth's interior can be classified into four parts: the solid inner core, the liquid outer 
core, the mantle, and the crust. We assume that the solid inner core is a homogeneous 
sphere (with density pi C ~ 1.29 x 10 4 kg m -3 and radius rj C ~ 1.22 x 10 6 m 22| ) and it 
is enclosed in the spherical liquid outer core (with the density of the fluid at the inner core 



boundary p oc ~ 1.22 x 10 4 kg m -3 22J). We do not consider any deformations (such as 
tidal and rotational deformations). We assume that the mantle and the crust are spherical 
shells with uniform densities, and their centers of figures are coincident; their gravitational 
influence on the inner core is negligible due to Newton's shell theorem. We do not consider 
any electromagnetic effects. Also, we assume that the friction (the friction coefficient 7) 
between the inner core and the outer core is proportional to the velocity of the inner core. 

By applying the generalized equation of geodesic deviation (Eq. (25) of [21I0, which 
includes the dilaton corrections to the standard equation of geodesic deviation [231], the 
translational motion of the inner core, relative to the rest part of the Earth, can approxi- 
mately be given by: 

~ _ 7 /i _ w 2f _ Aq{L k d k + 1)<9>, (1) 

where i and k indicate the space-time components of the parameters, Aq is the difference 
in dilatonic charge between the inner core and the rest part of the Earth (for the definition, 
see Eq.(5) below.), L is the rest separation between the Sun and the Earth (see Fig. CD), 
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Earth 




FIG. 1: A schematic cross section of the assumed configuration (not drawn to scale). As a dila- 
tonic wave passes along the Sun and the Earth (separated by L ~ 1.5 x 10 11 m), the inner core 
oscillates along the Sun-Earth line (amplitude /). The surface-gravity changes due to the inner core 
oscillations are searched for by superconducting gravimeters. We do not consider the influence of 
the Earth's rotation on the inner core for simplicity (see text). 



and is the dilaton field. The gravitational stiffness 24J and the friction coefficient are, 
respectively, given by 

^ ~ %*G Pic ~ Poc Poc « 1.8 x 10-V 2 « {2tt(4.1 h)- 1 } 2 , (2) 



3 

7 = « 2.3 x lO- 16 ^ b- 1 , (3) 



where G = 6.67 x 10~ n N m 2 kg -2 is the gravitational constant, mi c ~ 9.8 x 10 22 kg is the 
mass of the inner core, and r\ is the effective viscosity of the outer core. The stiffness due to 
the Sun's tidal force and the Moon's tidal force, which act to enlarge any displacement of 
the inner core from the center of the Earth, is negligible. 

We ignore the influence of the Earth's rotation on the inner core. The drag associated 
with the Earth's rotation takes its maximum value when it is balanced with the Coriolis 
force (~ 2rrii C ujiUJolo, where lor w 27r(24 h) _1 is the angular frequency of the Earth's rotation 
and Iq is a nominal magnitude of the inner-core displacement), which is about 30 % of the 



gravitational restoring force (m ic ci; 2 /o); for the liquid outer core, the drag would be less than 
the maximum value. The central force should be less than 2 % of the gravitational restoring 
force. 

The last term of Eq. ([T]) is due to the direct coupling to the gradients of the dilaton 
field 0. Aq is the difference in dilatonic charge between the inner core and the rest part 
of the Earth. For ordinary matter, the dilatonic charge can be obtained by summing over 



all the components n of the matter 



— I' 1 = ^; m Z q " ' Q * s defined as the relative strength 



of scalar to gravitational forces. For a body (mass M) composed by B barions with mass 
mj and (dimensionless) fundamental charge q b , and Z (~ B) electrons with mass m e and 
(dimensionless) fundamental charge q e , we obtain (Eq. (22) of [2l|) 

Bm b q b _ B 



M p 

M 



q b , (4) 



where /i = ^- is the mass of the body in units of baryonic masses. Therefore, Aq can be 



A^(-) -(-) }q b = A(^)q b , (5) 



given by 

M. 

where ic and rpe denote the inner core and the rest part of the Earth, respectively. The 
magnitude of A is about 10~ 3 or less (see e.g., Table 2.1 of pi). The upper limit on 
ql is determined to ~ 1.5 x 10~ 10 by tests of the equivalence principle [9(. Therefore, we 
obtain Aq < 1 x 1(T 8 . 

When a dilatonic wave propagates along the Sun-Earth line, the spectrum of the dis- 
placement along the line is given by 



at the resonant frequency: f = u /2n w 6.8 x 10~ 5 Hz, where Sh(f) is the power spectrum 
of the strain due to massless dilatonic waves: h = 20/ c 2 ; k is the wave number: k = u>o/c; 
and c is the speed of light. 

The intensity of a stochastic background of massless dilaton can be characterized by 



10|: Q DW (f) = (l/p c )(dp DW /d\ogf), where p c 



the dimensionless energy density 
3HqC 2 (8ttG) is the present value of the critical energy density for closing the Universe; Hq 
is the Hubble expansion rate: H = 3.2 x 10~ 18 h W0 s _1 ; h wo , ranging from ~ 0.5 to 1, is a 
dimensionless factor to account for different values of Hn. 
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The power spectrum of the strain is related to the dimensionless energy density by the 
following expression 

*W(/) = t£? S h(f) (7) 



3H, 



o 



for massless dilaton. Here the stochastic dilatonic background is assumed to be isotropic, 
unpolarized, and stationary lsf. 



IV. ESTIMATION OF THE UPPER LIMITS 



Coriolis acceleration splits the inner core oscillation into a triplet of periods (the Slichter 
triplet 25[|). To determine physical properties of the core, the three translational mode sig- 
nals (prograde equatorial, axial, retrograde equatorial) have been searched for in geophysics 



(.e.g. 



26 



27 



28 



29 



30 



31 



32 



33 Q, but the detection has not been confirmed yet. We use 



3 if ] to estimate the upper limit on the magnitude of the 



a recent analysis by Rosat et al. 
inner-core displacement. 

Rosat et al. examined gravity data obtained over the year 2001 from global observatories 
(in Canada, Australia, Japan, France, and South Africa) by applying a multistation stacking 
method [3jJ. They found no significant peaks that could originate from the translational 
motions of the inner core. The average noise level and the standard intervals were (4.2 ± 1.6) 
x 10 -12 m s~ 2 Hz -1 / 2 at ~ 7 x 10 -5 Hz in their analysis (Fig. 12(b) of [31|). This indicates 
that the upper limit on the amplitude of each translational mode is Sg = 5.8 x 10~ 12 m s -2 
Hz~ 1//2 at the frequency. This upper limit approximately corresponds to a displacement of 



the inner core: 



x 



2Gm, 



-5g « 0.11 mm Hz 1/2 , 



(8) 



where r P = 6371 km 



22] is the mean radius of the Earth. 



We assume that x represents the upper limit on the typical magnitude of the displacement 
With the value of x and Eq. (jHJ), we obtain the upper limit on the strain spectrum: 

27 



:V^)~l-5xl0- 20 ?7 Hz" 1 / 2 . 



Aq ■ c^L 2 uol/c 2 + 1 
From Eq. (J7D and the upper limit on the strain spectrum, we obtain 



(9) 



tt DW {fo)h 2 100 < 2.2 x 10-V- 



(10) 
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The effective viscosity r] is not well determined and estimates from various methods vary 
from ~ 10~ 3 Pa s to ~ 10 12 Pa s [34]. The Reynolds number (= p oc vri c /r], where v is the 
velocity of the inner core) is less than unity when 77 is larger than ~ 7 x 10 2 Pa s and 
the amplitude of oscillations x (= v/uq) is 0.11 mm. We obtain £2dw(/o)^ioo < 1 x 10~ n 
when the Reynolds number is unity and £|dw (/o)^ioo < 2 x 10 7 in the worst case when 
the effective viscosity is as large as 10 12 Pa s. When we use a recent estimate of an upper 
bound from nutation data [3j3|, r] < 1.7 x 10 5 Pa s, we obtain &Dw(fo)hi 00 < 6 x 10~ 7 . The 
expected upper limits for those different values of r\ are shown in Fig. [2j 



V. DISCUSSION 



By assuming the simple Earth model and using available geophysical results, we have 
estimated the upper limits on the dimensionless energy density of the stochastic background 
of massless dilaton. As one can see in the relation (ITU1) . the upper limits largely depend 
on the magnitude of the effective viscosity of the outer core, which is poorly constrained 
by experiments. Identification of the magnitude of the friction coefficient is crucial to this 
method. 

There are a number of active researches going on to determine the effective viscosity and 



other physical parameters of the Earth's interior. Some o 



them are laboratory measure- 



ments of the viscosity at high pressures and temperatures [36| , studies for neutrino oscillation 
tomography of the Earth's interior {3?]], and coincidence measurements with a laser strain- 



meter system and a superconducting gravimeter at the Kamioka Observatory in Japan [38]. 
With these researches employing new technologies, one can expect that our knowledge on 
the effective viscosity and dynamics of Earth's interior will be improved significantly in the 
near future. 

The dimensionless energy density of a massless background is constrained to fi/i 2 00 ^ 
10 -5 , in the frequency range of our interest, by the nucleosynthesis and resent measurements 
of the cosmic microwave background 39]. From the relation (flOj) . one can see that, if the 
effective viscosity is smaller than about 7 x 10 5 Pa s, the current sensitivity of this method 
is sufficient to reach the limit imposed by the astrophysical observations. As for the relic 
background of light dilatons, mentioned in the introduction, the intensity of the background 
is not constrained by the astrophysical observations and could compensate the weakness 



8 



10 





fammmmmmmmm. 














77 = 1 F 


a s 

























































SB SB 




























































m 




































o 




TU s Pas"" 
























































X 


•77— 700-F 


a-s 




iii 


ii 


III 


ii 


iii 


ii 



0.67 0.675 0.68 0.685 0.69 0.695 0.7 

x 10 

Freq uency[Hz] 



FIG. 2: Expected upper limits on ^Dty^-ioo f° r different values of 77. Asterisks, circles and diamonds 
indicate expected upper limits for r\ = 700 Pas (when the Reynolds number is unity), 1.7xl0 5 Pas 



(the upper bound estimated from nutation data 



35f|). and 10 Pas (the largest estimate given in 



34|), respectively (see text). 



of the dilaton coupling [101 ] . The possibility of testing the cosmological scenario has to be 
discussed in the future. 

We have assumed that the dilatonic waves are isotropic. However, if they are not isotropic, 
the magnitude of the inner-core displacement would depend on the location of the obser- 
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vatories. Such location dependence could be investigated through the Global Geodynamics 
Project network (GGP |40|) of superconducting gravimeters. 

Though this method is currently limited by the uncertainty in the Earth model, it is in 
principle sensitive to any composition-dependent waves at low frequencies. Other possible 
sources of such waves and their influence on the inner core have to be investigated. 

Geophysical sources of the excitation of translational motions of the inner core are not 



well known. Possible sources of the excitation are, for example, large earthquakes [41(, some 
dynamic flows 42J], and magnetohydrodynamic processes in the core 43]. The effects of 
earthquakes are thought to be very small [41] . If the translational motions are experimen- 
tally confirmed, all possible sources of the excitation have to be considered carefully before 
estimating the amplitude of the dilatonic waves. 

The sensitivity could be improved, for instance, by extracting correlated signals from the 
GGP network and applying an advanced data analysis method 44j . Also, the measurement 
sensitivity could be improved by further studies. 



VI. CONCLUSIONS 



We have proposed a new method of searching for dilatonic waves. In this method, Earth's 
surface-gravity changes due to translational motions of the inner core, excited by dilatonic 
waves, are searched for by using superconducting gravimeters. The main merits of this 
method may be as follows: (1) unlike the previous proposals, which intend to use gravita- 
tional wave detectors, this method employs the geophysical approach with superconducting 
gravimeters, (2) the sensitive frequency range is low (~ 7 x 10~ 5 Hz), in comparison with 
the previous proposals (~ 10 to 1000 Hz), (3) it is devoted to searching for the direct cou- 
pling of dilatonic waves to matter, and (4) it has a potential to search for anisotropy in 
composition-dependent waves, by observing anisotropy in the oscillation direction of the in- 
ner core through the GGP network. The major obstacle of this method is its dependence on 
the Earth model. This obstacle would diminish in time with the progress on understanding 
of the Earth's interior. 

Assuming the simple Earth model and configuration, we have obtained the preliminary 
upper limits on the stochastic massless dilatonic background at the frequency of fo ~ 7 
x 10~ 5 Hz. In the worst case when the effective viscosity (rj) is as large as 10 12 Pa s, 
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the upper limit on the dimensionless energy density (£lcw(/o)^ioo) is 2 x 10 7 . When the 
Reynolds number is about unity (77 ~ 7 x 10 2 Pa s), the upper limit is 1 x 10~ n . With 
the recent estimate of the upper bound from nutation data, 77 < 1.7 x 1CT 5 , we obtain 

^W(/o)^oo < 6x 10- 7 . 

This method is in principle sensitive to any composition-dependent waves at low frequen- 
cies. Other possible sources of such waves have to be investigated. We have focused on 
the massless dilaton. The possibility of searching for the relic background of light dilatons, 
predicted as a scenario in string cosmology, has to be discussed in the future. 

The sensitivity of this method could be increased with further studies, such as a develop- 
ment of coincidence measurements through the GGP network. With the increased sensitivity 
and refined Earth model, we could test predictions of unified theories and cosmological mod- 
els in a new window. 
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